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The results of a theoretical investigation on the low- velocity stopping power of the ions moving in a magne- 
tized collisional plasma are presented. The stopping power for an ion is calculated employing linear response 
theory using the dielectric function approach. The collisions, which leads to a damping of the excitations in the 
plasma, is taken into account through a number-conserving relaxation time approximation in the linear response 
O I- function. In order to highlight the effects of collisions and magnetic field we present a comparison of our analyt- 

ical and numerical results obtained for a nonzero damping or magnetic field with those for a vanishing damping 
or magnetic field. It is shown that the collisions remove the anomalous friction obtained previously [Nersisyan 
et ah, Phys. Rev. E 61, 7022 (2000)] for the collisionless magnetized plasmas at low ion velocities. One of 
major objectives of this study is to compare and contrast our theoretical results with those obtained through a 
novel diffusion formulation based on Dufty-Berkovsky relation evaluated in magnetized one-component plasma 
p£h \ models framed on target ions and electrons. 
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I. INTRODUCTION 

O 

' 55 ■ The energy loss of ion beams and the related processes in magnetized plasmas are important in many areas of physics such 
£>y as transport, heating and magnetic confinement of thermonuclear plasmas. The range of the related topics includes ultracold 
,-£h ' plasmas |Q]], cold electron setups used for ion beam cooling J2-01, as well as many very dense systems involved in magnetized 
target fusions lU, or inertial confinement fusion. This latter thermonuclear scheme presently advocates a highly regarded fast 
ignition scenario [6], based on femtolaser produced proton or heavier ion beams impinging a precompressed capsule containing 
7-H ■ a thermonuclear fuel 01 m it- Then, the magnetic field B values up to 10 10 G may be reached in the laboratory @]. Such a topic 
£> 1 is also of intense astrophysical concern [9]. These interaction geometries highlight low ion velocity slowing down (LIVSD) 
as playing a fundamental role in asserting the confining capabilities and thermonuclear burn efficiency in dense and strongly 
magnetized media. 

For a theoretical description of the energy loss of ions in a plasma, there exist two standard approaches. The dielectric linear 
response (LR) treatment considers the ion as a perturbation of the target plasma and the stopping is caused by the polarization 
of the surrounding medium [10-15]. It is generally valid if the ion couples weakly to the target. Alternatively, the stop ping 
is calculated as the result of the energy transfers in successive binary collisions (BCs) between the ion and the electrons II 16 - 
[l9ll . Here it is essential to consider appropriate approximations for the screening of the Coulomb potential by the plasma Oil- 
However, significant gaps between these approaches involve the crucial ion stopping along magnetic field B and perpendicular 
to it. In particular, at high B values, the BC predicts a vanishingly parallel energy loss, which remains at variance with the 
nonzero LR one. Also challenging BC-LR discrepancies persist in the transverse direction, especially for vanishingly small 
ion projectile velocity v when the friction coefficient contains an anomalous term diverging logarithmically at v — > ifTH 
[l2tl . In general when v is smaller than target electron thermal velocity v e , the ratio S(v)/v, where S(v) = —dE/dx is the 
stopping power (SP), usually monitors a linear stopping profile for highly ionized plasma with B = Il20ll or B ^ |Ht]. An 
alternative approach, particularly in the absence of any relevant experimental data, is to test various theoretical methods against 
comprehensive numerical simulations ll3l [l6l - tl8ll . The latter exhibit high a level of numerical noise at large magnetic fields, and 
in the v — > limit, while keeping a plasma coupling below unity, which is precisely the domain of many important applications 
of current interest. 

With this background we report on a theoretical study of energy loss of a slow-velocity ion in a magnetized classical plasma 
through a linear response approach which is constructed such that it conserves particle number. Broadly speaking, there are 
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two objectives of this paper. The first objective is to use the Bhatnagar-Krook-Gross (BGK) approach based on the Boltzmann- 
Poisson equations for a collisional and magnetized classical plasma which is treated as a one-component plasma (OCP). We use 
this approach to derive the dielectric function in a number-conserving manner and use this dielectric function to calculate various 
aspects of energy loss of an ion moving in the plasma. This is done in Sec. [II] We would like to see how number-conservation 
and damping (due to collisions) affect the stopping power of an ion in a low velocity range, i.e., for a slow ion. We should 
mention that for a colliisonal quantum plasma, e.g., a degenerate electron gas (DEG) without magnetic field, Mermin lf2lll and 
Das i22ll considered the equation of motion for a single-particle density matrix and derived the dielectric function in a number- 
conserving approach and in random phase approximation (RPA). The main part of our first objective is to see to what extent the 
BGK approach can address the SP of a slow ion in a magnetized one-component classical plasma. We will show that number 
conservation and collisions in such OCP have interesting and experimentally observable effects on low-velocity SP. Now, one 
may expect various collision mechanisms in a magnetized plasma. Because of their importance we have separately dealt with 
them in Sec. [Ill] In Sec. [IV] we present detailed calculations of low-velocity SP in the BGK approach. Although the LR is 
normally used to calculate the SP of a fast ion, we show that for the BGK model one can obtain useful and insightful results for 
low-velocity within the LR theory. 

A second objective of our paper is to compare and contrast our theoretical results with those obtained through a different 
method. The latter is specifically aimed at a low-velocity SP which is expressed in terms of velocity-velocity correlation and 
hence to a diffusion coefficient. We refer to Dufty and Berkovsky (DB) [23] for an exposition of this method. The later part 
of Sec. [IV]contains a brief account of the DB relation. Marchetti et al. [24] and Cohen and Suttorp l25ll have calculated the 
relevant diffusion coefficients in a magnetized plasma through a hydrodynamical and kinetic approaches, respectively. Recently 
Deutsch et al. in Ref. J26ll have suggested an alternative approach for a calculation of low-velocity SP via the DB relation and 
by employing the diffusion coefficients for a magnetized OCP 1I24II25I1 . featuring either the target electrons or target ions. Since 
our main theoretical results are obtained in a kinetic equation approach we have decided to devote Sec.[V]to an appraisal of these 
two approaches. In that section we discuss results for an unmagnetized and magnetized plasma in the two approaches. This 
serves to highlight the merits of the kinetic approach versus those in a hydrodynamical approach. 

Section[VT]contains some discussion and outlook. Appendix[A]can be consulted for an integral representation of the dielectric 
function. The Coulomb logarithm where the dynamical polarization effects are neglected is briefly discussed in AppendixlBl 



II. LINEAR RESPONSE FORMULATION 



In this section we consider the main aspects of the linear response (LR) theory for the ion-plasma interaction in the presence 
of an external magnetic field. Within the LR, the electron plasma is described as a continuous, polarizable medium, which is 
represented by the distribution function of the electrons /(r, v, t). The evolution of /(r, v, t) is determined by the kinetic and 
Poisson equations. Usually only a mean-field interaction between the electrons is considered and hard collisions are neglected. 
This is valid for weakly coupled plasmas where the number of electrons in the Debye sphere No = 47rno e Ag = 1/e > 1 is 
very large. Here e is the plasma parameter, n$ e and A e = (/csr e /47rno e e 2 ) 1 / 2 are the equilibrium density and the Debye length 
of electrons, respectively. 

We consider a nonrelativistic projectile ion with charge Ze and with a velocity v, which moves in a magnetized collisional 
and classical plasma at an angle d with respect to the constant magnetic field B. We ignore any role of the electron spin or 
magnetic moment due to the nonrelativistic motion of the ion and the plasma electrons. We shall consider here the limit of heavy 
ions and neglect recoil effects. The strength of the coupling between the moving ion and the electron plasma is given by the 
coupling parameter 



Z = 



Z/N D 



(1) 



Here v e is the thermal velocity of the electrons. The derivation of Eq. (fl~|l is discussed in detail in Ref. H27I1 . The parameter Z 
characterizes the ion-target coupling, where Z -C 1 corresponds to weak, almost linear coupling and Z > 1 to strong, nonlinear 
coupling. 

Let us now specify the kinetic equation for the collision-inclusive classical magnetized plasma. The effect of collisions on 
the dielectric properties of the plasma is included, in a number-conserving approximation, through a relaxation time r = 1 /j, 
where 7 is the collision frequency [28]. For r — > 00 the collision-inclusive kinetic equation reduces to the collisionless Vlasov 
equation. Thus we consider the kinetic equation of the collisional plasma within relaxation-time approximation (RTA) in which 
the collision term is of the Bhatnagar-Gross-Krook(BGK)-type [28], 
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where the collision frequency 7 is a measure of damping of excitations in the plasma, and 



n e (r, t) = J f (r, v, t) dv, n 0e = J 



fo (v) dv. 



(3) 



Here n e (r, t) is the density of the electrons, /o(v) is the equilibrium distribution function of the electrons in an unperturbed state. 
For instance, the distribution function /o(v) of the plasma electrons may be given by the Maxwell distribution. The right hand 
side of Eq. (O is the collision term in a relaxation-time approximation which was introduced by BGK in a number-conserving 
scheme. It is easy to see that this form of collision term conserves the total number of particles. E is a self-consistent electric 
field (see below) and B is treated as an external magnetic field, t — I/7 is the relaxation time. 

For a sufficiently small perturbations (LR treatment) we assume / = f + fx, n e — n 0e + n le , with 



nie (r, t) = J fi (r, v, t) dv, 



and the linearized kinetic equation becomes 

dfi , dh 
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(5) 



Here b = e z — H/B is the unit vector along the magnetic field, 51 e = eB/m e c is the cyclotron frequency of the electrons, 
E = — Vy, if is the self-consistent electrostatic potential which is determined by the Poisson equation 



V 2 ip = -471-po (r, t) + 47re J f x (r, v, t) dv, 



(6) 



where po is the density of the external charge. 

We solve the system of equations (|5j and (|6]l by space-time Fourier transforms. Because of the cylindrical symmetry of 
the problem around the magnetic field direction b, we introduce cylindrical coordinates for the velocity v = e x v± cos 9 + 
e y v± sin6* + e z W|| and the wave vector k = e x k± cosi/}+e y kx sin %/) + e z ku, where the symbols || and _L denote the components 
of the vectors parallel or perpendicular to the external magnetic field, respectively. We next introduce the Fourier transforms of 
/i(r, v, t), nie(r, t) and tp(r, t) with respect to variables r and t, /ik,w(v), «i e (k, uS) and <ys(k, uj). Then the linearized kinetic 
equation ((5]) for the distribution function becomes 



dfikojjv) 
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Assuming axially symmetric unperturbed distribution function, /o( v ) = /o(|^|| 
solution is 



Eq. (0 can be formally integrated and the 
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d9', 



where the lower limit of 9' -integration is chosen so as to take the integrand vanish. Here 

U (9 r ) = (k\\v\\ -cj-i-y) {9' -9) + k±v± [sin (9' ~ tp) - sin (0 - tp)} 



(9) 



The 9' -integration in Eq. (JHJ can be performed using the Fourier series representation of the exponential function i29ll . After 
straightforward integration we obtain 



fi-kw (v) = -e" 



^n(9^)j n{ze) 



hnie (k, ui) + — ip (k, u) 

n 0e m e 
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where z e = k±v±/il e , J n is the Bessel function of the nth order. It should be emphasized that Eq. (fTOb is a formal solution of 
the linearized kinetic equation because the Fourier transformed electronic density n le (k, uS) remains unknown. We now perform 
v-integration in Eq. ( fTOb and solve the obtained algebraic equation with respect to the quantity ni e (k, uj). Substituting this 
quantity into Fourier transformed Poisson equation finally yields 

47rpo (k, uj) 

Pfcu =75 j- r, (11) 

k 2 e M (k,w,7) 

where e^f (k, u, 7) is the collision-inclusive longitudinal dielectric function of the magnetized plasma which is given by 

e M (K^) = l+ ( " + i l n f' UJ ' l) : 1] (12) 
UJ + 17Q (k,w,7) 

with e(\s.,oj,j) = e e (k,o; + ij), Q(\l,uj,j) = Q e (k,u; + ij) and 

uj 2 2ir f°° f°° 

e e (k, w) = 1 - -f / dv|j / u_i_<fo_i_ (13) 
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(14) 



Here e e (k, cj) is the usual longitudinal dielectric function of the magnetized collisionless and purely electron plasma (see, e.g., 
Ref. l28ll ) and uj e — (47rno e e 2 /m e ) 1 / 2 is the plasma frequency of the electrons. Similarly, Q e (k, uj) refers to the electron plasma. 
The dielectric function Em (k, w, 7) given by Eqs. (fT2l-(fT4l has been obtained in the BGK approach which is number-conserving. 
Note the exact relation Q e (k, 0) = 1 which holds independently of the initial distribution fo if the latter is normalized to the 
unperturbed electronic density no e , see the second relation in Eq. (01. 

It is well known 1 28] that the usual relaxation-time approximation can be obtained from Eq. (O if the collision term is written 
as — l{f ~ fo) and is equivalent to replacing uj by uj + 17 in the collisionless dielectric function e e (k, uj). This procedure is 
inadequate because it does not conserve the local particle number and does not lead to the Drude behavior at long wavelengths 
(k — s- 0). This is remedied in the BGK approach. The k —> case of the number-conserving dielectric function is also of 
interest. Noting that Q e (0, uj) — 0, from Eq. (fT~2T > we find 



r„U') = p^±M + ph|M (15) 



with 



s ± (uj) = 1 + ^ ( " + Z7) a , £|| H = 1-— ^— . (16) 



The above results are a generalization of the Drude dielectric function for magnetized plasmas. Equations (TT3T > and ([Tol l are 
known also as a "cold" plasma approximation (see, e.g., Ref. 12811 ) and can be alternatively obtained from Eqs. (IT2~b-(fT~4T> assum- 
ing initial distribution function /o(v) = no e S(v). For a simplicity we shall call the expressions ( fl~3T > and dT~4l > as a Bessel-function 
representation of the dielectric function. For many practical applications, however, it is important to represent the dielectric func- 
tion in an alternative but equivalent integral form, see Appendix|A]for details. 

Let us now specify the initial distribution function fo of the electrons. We consider the Maxwell isotropic distribution function 

t ( \ n °e I V l +v2 A n-n 

fo ^ v±)= ^^ exp {-^r)> (17) 

where v e — (fe^Te/me) 1 / 2 is the thermal velocity of the electrons. The collision-inclusive dielectric function then reads 

e (k, W) 7) = 1 + ^ t F i ( k > w ) + iF 2 (k, uj)] . (18) 



Here 



Fi(k,w) = l+ V ^ n A n (fie) [uQ (x n , y) - 7-F {x n , y)} , (19) 
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FIG. 1: (Color online) Generalized Drude (solid line), collisionless (dashed line) and number-conserving BGK (dotted line) energy loss 
functions vs uj/ui e for Q e = uj e , 7 = 0.1u e and ku A e = k±X e — 0.1 (left panel), ku A e = k±X e = 0.4 (right panel). 



are the real and imaginary parts of the generalized dispersion function of the collisional magnetized plasma, respectively, and 
x n = (uj + nf2 e )/|fc|| \v e , y — j/\kn \v e , f3 e = k\a 2 , a e — v e /£l e is the cyclotron radius of the electrons, A„(z) = e~ z I n (z), 
I n (z) is the modified Bessel function of the nth order. Here we have introduced the generalized Fried-Conte dispersion functions 
for the collisional plasma 

G{x ' y) = 7^L {t-xf + t? ' (21) 



F(x,y) = ^L / . (22) 



,2 



27T J_oo (t — x) +y 

At vanishing 7 (at y — > 0) these functions become the usual Fried-Conte dispersion functions [30] of the collisionless plasma 

'W-lsL—' ,23) 

n -x 2 /2 



where 



F{x,Q) = ^-xe- x > 2 . (24) 

The function Q(k, £^,7) which determines the dielectric function (fT2l is evaluated by inserting Eq. ( fTTI i into Eq. ( fl~4l . It is 
easy to see that for the Maxwell isotropic distribution function (fTTT i the quantity Q(k, ui, 7) is given by 

Q (k, w, 7) = £e (k ',t n^r 1 = Fl ( k ' w ) + lF2 ^ ^ ' (25) 
e e (k, 0) - 1 

e e (k, 0) = 1 + ^ 26 > 

is the static dielectric function which is not affected by the external magnetic field. 

For ion stopping considerations, it is worth defining the energy loss function (ELF) Im[— l/e(k, uj,j)]. Figure [T] shows 
Drude, collisionless, and BGK energy loss functions vs scaled frequency Ld/uj e when fci|A e = k±\ e = 0.1 (left panel) and 
fen A e = k±X e = 0.4 (right panel) for £l e = uj e , 7 = 0.1u) e . As has been mentioned above at small momentum k — >• the 
BGK energy loss function reproduces the Drude energy loss function. And this is seen on the left panel of Fig. [T] Also at long 
wavelengths (i.e., at small k) the BGK energy loss function is broadened due to the damping compared to the ELF with vanishing 
damping. 

The collision-inclusive dielectric function allows both physical insight and useful numerical estimates of the influence of the 
collisions on energy loss. In an unmagnetized and degenerate electron gas the predicted effect is a shorter lifetime and smaller 
mean free path of the plasmons resulting in considerable modifications of the ELF 113 1143511 . For the stopping of a single ion, 
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the broadening of the plasmon peak with increasing 7 shifts the threshold for the energy loss by plasmon excitation towards 
lower projectile velocities. This increases the SP at low projectile velocities, compared to the collisionless result 13 1143411 . The 
situation with a present case of a classical and magnetized plasma including the collisions may be quite different although 
collisional broadening of the ELF occurs also in this case. This situation will be further discussed in Sec. lIVI 



III. COLLISION FREQUENCY IN A MAGNETIZED PLASMA 

In Sec.UJthe effect of the collisions in a magnetized classical plasma has been introduced in the dielectric function through 
a phenomenological but number-conserving collision term within the LR theory. The model collision frequency 7 in solids and 
plasmas can be determined experimentally or, alternatively can be calculated theoretically. For instance, in some investigations 
of ion stopping in solid targets in the absence of a magnetic field, 7 was determined by fitting — Im[£ (0, w, 7)] to experimental 
optical energy loss functions (see, e.g., Refs. IT331 [34(1 and references therein). In addition the model relaxation time r = 1/7 
can be estimated from the experimental data of the dc conductivity or the mobility in a plasma either with or without external 
magnetic field. It should be emphasized that in general there are a number of physical mechanisms which may contribute to the 
damping parameter 7. And contribution of each mechanism depends strongly on the specific plasma conditions. We have not 
attempted here to evaluate the damping parameter from first principles in the most general case but regard it rather as a model 
parameter. In principle 7 can be calculated to varying degrees of approximations which may allow us to see how the SP depends 
on the target properties and the magnetic field through their influence on 7. 

In this section we briefly consider a fully ionized and a weakly coupled plasma where the contributions of the Coulomb 
collisions to the frequency 7 may play a dominant role. This frequency, in our case, is determined by electron-electron (e-e, 
7ee) and electron-ion (e-i, -f e i) Coulomb collisions (if we do not consider impurities). Thus, in contrast to Sec.HU we deal with 
a two-component electron-ion plasma (TCP) accounting for the dynamics of plasma ions. The total effective frequency, in the 
limit of a weakly coupled plasma, can be approximated as a sum of e-e and e-i collisions, 7 = 7 ee + 7^. In the absence of a 
magnetic field the theory of Coulomb collisions in a plasma has been formulated by Spitzer |36ll (see also |H3]). In the last four 
decades or so the theory has been further developed and extended. The recent book 1 38] summarizes the results obtained during 
last four decades. However, to our knowledge, the relaxation processes in a magnetized plasma have not been studied in as much 
detail as in an unmagnetized plasma, and only several theoretical attempts exist for this case ll39l - l46ll (see also the references 
therein). For a classical plasma more complete expressions for the collision frequencies valid at arbitrary (but non-quantizing) 
magnetic fields have been derived by Ichimaru et al. and Matsuda 141H44I1 . and by Montgomery et al. 14511 and by Silin H4al 
with and without allowing for dynamical polarization effects in plasma, respectively. 

In Refs. lf4lL l42tl only e-i relaxation is considered. The generalization to the e-e case is straightforward. The final result is 
summarized by a formula 



8V27n^e 4 n a ?7 ea 

lea = — 5 ; In A eQ , (27) 

3m e m a v^ a 

where a = e, i indicates the plasma species, q e = — 1, = Z^ Z^e is the charge of plasma ion, rj ei = 1, rj ee = 2 1 / 2 , 

1 f dk r- Ge (^)G a (^)f d ^ (28) 



2(2w) 3/2 v 3 ea J J -00 fc2 fc4 | eei(k)W) | 



G a (k,w) 



A„ (k 2 ± a 2 a ) exp 



(lu — nf2 Q ) 



(29) 



v 2 a = k B T a /m a , v\ a = v\ + v%, v~l = v~ 2 + v~ 2 , a a = v a /Q a , A n (z) = e~ z I n (z). Also T if rrii, Vi and Q t = 
ZieB/rriiC are the temperature, the mass, the thermal velocity and the cyclotron frequency of the plasma ions (aj is the cyclotron 
radius), respectively. The quantity In A ea is the generalized Coulomb logarithm for a magnetized plasma. Here £ e i(k, ui) is the 
longitudinal dielectric function of a magnetized and collisionless electron-ion TCP (see, e.g., Ref. JH). The limit of the 
vanishing magnetic field in Eq. ( 1291 is not trivial. An alternative but equivalent integral form for the function d29| > allowing 
easily the limit of the field-free case is derived in AppendixlAl see Eqs. (IA14t and ( |A151 l. In Eq. d28| i the dynamical polarization 
effects are included in a dielectric function £ e i(k, oS). These effects guarantee the convergence of the k-integration in Eq. (|28l i 
at large distances or at small k. But an upper cutoff fc max = 1/ r m - m (where r m ; n is the effective minimum impact parameter) 
must be introduced in Eq. ( |28] i to avoid the logarithmic divergence at large k. This divergence corresponds to the incapability 
of the linearized kinetic equation to treat close encounters between the plasma particles properly. Also it should be emphasized 
that for the e-i collisions there are two specific frequencies 7i e = Zi^ e i and 7 ej ; which describe the relaxation of ionic and 
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electronic temperatures to their equilibrium values, respectively l4ll l42tl . Thus the total e-i collision frequency is given by 

lei = lei + lie = {Zi + l)lei- 

To estimate the range of variation of the collision frequency with increasing magnetic field consider now some particular 
cases. At vanishing magnetic field from Eqs. ( lAT4b and ( |A15t we obtain G a (k, w) = \/k) exp(— u 2 /2k 2 v 2 t ). In this limit 
we denote -f ea = 70, ea with A ea = Ao, e a and from Eq. ( l28b we find 



2 r k ™* dk f°° e~ u2 ' 2 u 2 du 
7T Jo k J \ £o ei kv ea u)\* 



lnA , ea = \/~ I —I t-t^ -5 (30) 



the Coulomb logarithm in the absence of a magnetic field. Here £o, e i(&, is the usual longitudinal dielectric function of the 
electron-ion unmagnetized TCP without collisions. Equation ( f30b has been derived by Ramazashvili et al. IB7I1 . It involves 
the dynamical polarization effects through the dielectric function So,ei(k, ui) and requires only an upper cutoff fc max = 1/Vmin 
in a Fourier space. The Spitzer formula is recovered assuming £o, e i(k,Ld) = 1 in Eq. ( f30b and introducing a lower cutoff 
&min = 1/Ad- Then performing u-integration in Eq. ( f30b we obtain the usual Coulomb logarithm with Ao jeQ = Ad/V 6 



' ea,mm 



generalized for electron-ion plasmas (see, e.g., Ref. ll35lP . Here A^ 2 = A~ 2 + A~ 2 while r eQ m i n = max[ALea; Adb]> where 
A« = v a /uj a and uj a = (47r7i 0Q (j ,2 e 2 /TO C[ ) 1 / 2 are the Debye screening length and the plasma frequency for plasma species a, 
respectively. X^ ea denotes the usual Landau length X^ ee = Xhei/\Zi\ = e 2 /3/csT e and Adb = fi/2^/m e fcsT e , the electron de 
Broglie wavelength, taking care of the intrinsically quantum behavior of the high-temperature plasma in the short range limit. 

In the opposite case of an infinitely strong magnetic field (j ea — loo.ea with A ea = Aoo eQ ) from Eqs. dA14b and ( IA15b 
one obtains G a (k, u) = exp(— uj 2 /2k 2 v 2 t ). Inserting this formula into Eq. ( f28b it is straightforward to show that the collision 
frequency in a strong magnetic field is the half of the frequency io,ea- Thus 

In Aoo^a = - 111 A 0j ea- (31) 

The collision frequencies have been also investigated by some authors using the Fokker-Planck kinetic equation with the Lan- 
dau integral of the collisions l45ll46ll . As stated above this approach neglects the dynamical polarization effects. In AppendixlBl 
we show briefly that starting with Eqs. d27| i and d28| i where £ e i(k, uj) is set = 1, one arrives at the expressions derived by Silin 
ll46ll . To understand the importance of the dynamical polarization effects which are neglected in Eqs. (IB2I) and dB6b (and also 
in the formula of Spitzer) we note that the kinetic equation in the form of Landau accounts for only close (almost) Coulomb 
collisions, thus completely neglecting long-range wave-particle interactions. It has been shown previously that in the absence of 
magnetic field corrections A eQ to the standard Coulomb logarithm \n(X^/r ea , m i u ) arising due to these interactions may be of 
the same order as the leading term [37]. However, this effect is crucial only for e-i interactions and wave-particle interactions are 
not expected to make any essential change in the rate of e-e Coulomb collisions. For instance, in an anisothermic electron-ion 
plasma low-frequency ion-acoustic waves may provide an effective mechanism for electron-ion interactions which leads to an 
enhancement of the standard Coulomb logarithm. 

Similar long-range and low-frequency collective effects are responsible for a strong enhancement of the collision rates in a 
magnetized plasma. In this case it has been shown l4ll l42tl that e-i collision rate contains, in addition to the contribution from 
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TABLE I: The Coulomb logarithm J32b for electron-proton plasma normalized to In Ao, e i for some values of the scaled magnetic field B/B 3 
and for different values of the cutoff parameter £ = AD/V e i,mm = 10, 10 2 , 10 3 . 



B/B s 


0.1 


1.0 


10.0 


10 2 


10 3 


10 4 


10 5 


10 6 


( = 10 


0.05 


0.53 


3.11 


6.79 


10.55 


14.31 


18.07 


21.82 


£ = 10 2 


0.03 


0.27 


1.55 


3.40 


5.28 


7.16 


9.03 


10.91 


£ = 10 3 


0.02 


0.18 


1.04 


2.27 


3.52 


4.77 


6.02 


7.27 



close Coulomb collisions (IB61 . a term which in our notations is given by 

lnA*j = - In -M e z [(1 + z) K (z) - zK x [z)\ , (32) 

where z = l/£ 2 , C — A e /a e = is the scaled magnetic field, Kq(z) and K\{z) are the modified Bessel functions of the 

second kind. The Coulomb logarithm ( l32t depends essentially on the magnetic field. At large and small magnetic fields Eq. (l32b 
behaves as lnA* ; ~ ln(f2 e /w e ) and lnA^ ~ il e /ui e , respectively, and vanishes at il e — > 0. Also we note the large factor 
ln(mj/m e ) in Eq. d32l which diverges at mi — » oo and appears due to strong electron-ion interactions via collective plasma 
waves and is typically > 10. However, Eq. (IB6l l which accounts for only close Coulomb collisions does not contain such term 
and remains finite at rrii oo. 

The results of the numerical evaluation of Eqs. (IB21 i and ( IB6l l are shown in Fig. [2] This figure shows the ratio In A ea / In Ao jeQ 
as a function of the scaled magnetic field B/B s for e-e (left panel) and electron-proton (right panel) collisions and for different 
values of the cutoff parameters (£ e = AD/V e e,min ( e_e ) an d £ = AD/^ei.min (e-i)). The quantity B s introduced above is 
B s = mcv e /e\r>- Also for an electron-proton plasma the normalized Coulomb logarithm In K* ei j In A u>e j, Eq. d32l , is given in 
Table Ufor some values of the scaled magnetic field. It is seen that the magnetic field may essentially increase the collision rates 
in plasma compared to the field-free ones and this is more important for electron-ion collisions. As discussed above, Eqs. (IB2t 
and ( IB6b and hence the results shown in Fig. [2] account for only unscreened Coulomb collisions neglecting dynamic polarization 
effects. In a vanishing magnetic field these effects are important for e-i collisions and the situation with B ^ requires further 
investigations, in particular for e-e collisions. In the e-i case a major contribution is expected from low-frequency collective 
modes given by Eq. d32l . Table U shows that this contribution exceeds the Coulomb logarithm In A e j at B/B s ~ 10, 10 2 , 10 3 
depending on the cutoff parameter. Also it should be emphasized that the validity of the regime OTb of a classically strong 
magnetic field (the domains B/B s > 10 4 and B/B s > 10 6 in the left and right panels of Fig. [2l respectively) requires the 
condition Kl e < hsT? (or B < B c = (mc/ eh)kBT e ). Thus the results shown in Fig.|2]are valid up to B c /B s = ksT e /fujj e . 
Clearly the realization of the regime PIT ) requires high temperatures and low densities and the enhancement of the collision 
rates at B ~ 1 — 10 2 _B S may not be accessible under certain conditions. However the recent analysis shows ll39ll (see also the 
references therein) that in a quantizing magnetic field with B > B c the field-dependence of the collisional rates becomes even 
stronger and the enhancement of j ea shown in Fig. [2] may turn even more significant although the classical expressions d27l ), 
( |28| >. (IB2b and ( IB6I ) are invalid in this regime. 



IV. LOW- VELOCITY STOPPING POWER 

In this section, with the collision-inclusive dielectric function derived in Sec. [II] we consider the stopping power (SP) of a 
low-velocity ion moving in a magnetized plasma for an arbitrary angle with respect to the magnetic field. The regime of low 
velocities is of particular importance for some physical situations, e.g., for electron cooling processes [3] and for magnetized 
target plasma fusion researches [5]. Previously the SP in a magnetized plasma at small ion velocity has been investigated by 
employing linear response (LR) theory [11] and Dufty-Berkovsky relation 12611 . The latter approach (see below) reduces the 
problem to a determination of the diffusion coefficient of the magnetized plasma. In Ref . [11] it has been shown that in the 
presence of a magnetic field and in the absence of collisions, the friction coefficient contains an anomalous term which diverges 
at v — > like ln(v e /v) in addition to the usual constant one while the hydrodynamic approach of Ref. 12611 does not contain such 
term. We shall comment on this feature in this section. 

The stopping power S of an ion with charge Ze and velocity v is defined as the energy loss of the ion in a unit length due to 
interaction with a plasma. From Eq. ( flTT i it is straightforward to calculate the electric field E = — Vy> (or E(k, uS) = — zkyj(k, to) 
in terms of Fourier transforms), and the stopping force acting on the ion. Then, the stopping power of the projectile pointlike ion 
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becomes (see, e.g., Refs. [10-12]) 



S = 



2ir 2 v 



dk ^~ J im- 



-1 



(33) 



£m (k, k • v,7) 

For the friction coefficient we have to consider S, given by Eq. d3~3l in a low-velocity limit, and thus the dielectric function 
(IT2b with (TT~8T > and the functions -Fi(k, and ^(k, w) given by Eqs. ( lAlOt and (I A lib , when io = k • v. Now we have to 
write the Taylor expansion of Eq. ( fT2l for small ui = k ■ v. Using expressions ( fT2b . ( fTSI l, d25"l l. (|26*T i, ( IAIOI i. and ( IA1 II ) for the 
collision-inclusive dielectric function at to — > we obtain 



Im 



-1 



k\l 



e M (k,w, 7 ) (k 2 X 2 e + l) 2 v e 



-X(t)-ct 



(34) 



where <r = j/kv e = v jk\ e , v = 7/w e . The function is determined by Eq. JA9I ). It should be emphasized that Eq. d34"i > does 
not contain any logarithmic singularity at vanishing fen — > as for the case of collisionless magnetized plasma, see Ref. fiUl . 
This singularity which leads to an anomalous friction in a magnetized plasma has been removed here due to the collisions and 
the factor e~ qt in Eq. d34l ) guarantees the convergence of the t integration at kn — > 0. Thus from Eqs. ( f33l > and d34l i we obtain 
usual (linear with respect to v) friction law 



S(0) = 

where ft($) is the dimensionless friction coefficient, 

_ /nN f K k 3 dk 

n W = / 771 

Jo (k 2 + If 



i*k\1 v, 



ipx(k) cos 2 + -^(fc) sin 2 z? 



(35) 



(36) 



Here k = k max _X e an d $ is the angle between v and b. In Eq. d36b we have introduced a cutoff parameter fc max = l/r m i n 
(where r m ; n is the effective minimum impact parameter) in order to avoid the logarithmic divergence at large k. This divergence 
corresponds to the incapability of the linearized kinetic theory to treat close encounters between the projectile ion and the plasma 
electrons properly. For r m j n we thus use the effective minimum impact parameter of classical binary Coulomb collisions which 
at low-velocities of the ion reads r m i n = \Z\e 2 /mv 2 . It is seen that the parameter k = 4irno e Al/\Z\ = 1/Z ^> 1, where Z is 
determined by Eq. ([1) at v <C v e . Also the other quantities in Eq. ( T3"6j > are 



Ipn (k) 



exp 



2tf 



sin 2 (CO - 2vt 



*n (kQ(t)) 



dt 



tr (a) 



(37) 



with n = 1,2, C = A e /a e = n e /u e , Q(t) = y/2tT ((t), T 2 (t) = 1 - (sini/i) 2 , $i(x) = x- 2 $(x), $ 2 (x) = 2eii(x) - 
x~ 2 $(x). The function $(x) is determined by Eq. ( IB5b - 

In many experimental situations, the ions move in a plasma with random orientations of $ with respect to the magnetic field 
direction b. The friction coefficient appropriate to this situation may be obtained by carrying out a spherical average over $ of 
ft (tf) in Eq. (ED We find 



k 3 dk 

(k 2 + if 



[Vi (*) + ^2 (k)} 



(38) 



Let us analyze the general expression ( l36l l for some particular cases. For instance, at vanishing magnetic field (C — > 0) using 
the relation Q(t) ~ \fW&(,t 2 at C 0, one finds 



(39) 



where = e 2 erfc(z), erfc(z) is the complementary error function. In this case the friction coefficient is isotropic and 

becomes 



4 dfc 

k(2k 2 + iy 



(40) 



Here 1 /po — \J~2k. In addition at vanishing damping, i.e. at v — > 0, A(vk) — > 1 and we recover the usual low-velocity stopping 
power in an unmagnetized collisionless plasma with a friction coefficient (see, e.g., ITlll47ll ) 



fto(tf) = ±U(k) ee 1 



ln(l 



(41) 
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At strong magnetic field (£ —> oo), the plasma becomes highly anisotropic and the friction coefficient depends essentially on 
the angle d. For an evaluation of the functions ipi(k) and ip2(k) we note that Q(t) — > \/2t and T(£t) — > 1 as ( ^ oo. Then 
substituting these relations into Eq. d37] > and after integration by parts one obtains 



with a = v /y/2k, and 



B(z) 



^ (k) = lA(a)+a 2 B(a)-^=, 

yj 2 ( k ) = (l-a 2 )B(a)-±A(a) + -j= 



00 dt ... 2z 
A{t) = -= 
t Jtt 



In 



1 e" 



dt. 



(42) 
(43) 

(44) 



Then the friction coefficient at infinitely strong magnetic field reads 

1 f°° dk 



ftoo(tf) = 

+ cos 2 ■d 



2.L, fc(2fc 2 + l) 2 



l-u 2 k 2 )B(uk)--(A(uk)-^ 
2 V Jit 



(45) 



A (uk) -=vk + 2v 2 k 2 B (vk) 

J 7T 



Similarly for the angular averaged friction coefficient we obtain 



(ftooO?)) = 3 



pi. 



B(vk)dk 
k(2k 2 + if 



(46) 



The function B{z) involved in Eq. d45l > at small z behaves as B(z) ~ ln(l/z) — C/2, where C — 0.5772 is the Euler's 
constant, and diverges logarithmically when z — > 0. Using asymptotic behavior of this function it is straightforward to calculate 
from Eq. (|45T > the friction coefficient at vanishing 7. In this limit and in the leading order we obtain 



ftoo(0) = j < sin 



where is given by Eq. (HTl l. and 



2 t3 



In 



\/2a 



C+l 



£7 (k) + Ui (k) 



17 («) cos 2 1? 



(47) 



C/i (k) = f7 (k) In « - - [In 2 (k 2 + l) - 2 In (k 2 + l) 



1 l k- 
- -Li 2 
2 



(48) 



Here Li2(z) is the dilogarithm function. Note that at large n 3> 1, which is a requirement of a weak ion-plasma coupling, the 
functions Ui(k) and U(k) can be approximated by U\(n) ~ In 2 k — 7r 2 /12 and 17 (k) — 21nK — 1, respectively. It is seen 
that the first term in Eq. d47b diverges logarithmically at vanishing 7. It can be shown that the general expression d36l l with 
(f37T > for the friction coefficient derived for arbitrary but finite magnetic field behaves similarly. This is a consequence due to 
the magnetic field since the field-free result d40i > remains finite as 7 — > (see, e.g., Eq. d4Tl )). The divergent term in Eq. d47l ) 
vanishes, however, when the ion moves along the magnetic field (d = 0). Then the friction coefficient is solely given by the last 
term of Eq. d47i >. In addition, the friction coefficient Eq. d47b for strong magnetic fields shows an enhancement for ions moving 
transverse (i9 = 7r/2) to the magnetic field compared to the case of the longitudinal motion (d = 0). This effect is in agreement 
with particle-in-cell simulation results 

As stated in Introduction we shall now make contact with a different method. It has been shown by Dufty and Berkovsky llzll 
that the low-velocity SP of an ion in a plasma is related to the diffusion coefficient D through 



S(v) 



D 



(49) 



As in Ref. [26] we consider D to the self-diffusion coefficient in a magnetized classical one-component plasma. From Eqs. ( f35l >- 
d37l i we can relate the friction coefficient 7?.(i9) to the diffusion coefficient D through Eq. j49l . Cohen and Suttorp [25] have 
calculated parallel (to the magnetic field) diffusion coefficient Dm. These authors, like us but unlike Marchetti et al. [24], have 
used a kinetic equation method. At vanishing damping (7 — > 0), it can be shown that D11 obtained from Eqs. (f35T> - (T37T> and 
for d = coincides with the result of Cohen and Suttorp lf25Tl . In part icular, at 7 — > 0, it is found from Eqs. fiH and d47l ) 
that 7?-o(0)/7?. 00 (0) = -Doo,||/-Do,|| = 2/3 in agreement with Ref. [i25j[. Here -D00JI an d -Do,j| are tne parallel self-diffusion 
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FIG. 3: (Color online) The friction coefficient TZ(-d) vs the scaled magnetic field Q e /uj e (left panel) and damping parameter "//oj e (right panel) 
for i? = (solid line), $ = 7r/4 (dashed line) and i9 = n/2 (dotted line), k = 10, 'y/oJe = 0.1 (left panel), = ui e (right panel). The line 
with symbols corresponds to (7£(i?)). 



coefficients at infinite and vanishing magnetic fields, respectively. However at finite 7 and for 1} = 0, comparing Eq. d40b with 
d45l i we conclude that the simple relation cited above is not obeyed in general, due to damping. 

As an example, we show in Fig.[3]plots of the dimensionless friction coefficient given by Eq. J361 , as a function of 

the scaled magnetic field £l e /uj e for model parameter 7 (prfu) e — 0.1) (left panel). The right panel of Fig. [3] shows 7£(#) as a 
function of the scaled damping parameter 7/w e f° r = w e (i.e. for a given magnetic field). It is seen that the low- velocity SP 
increases with an increase in the angle -d and also with the magnetic field. In the latter case the SP asymptotically tends to the 
value given by expression d45l l. In the opposite limit of a weak magnetic field the friction coefficient tends to the value given 
by Eq. d40l which is independent of the angle d. Also Fig. [3] shows that the friction coefficient decreases with damping. It is, 
interestingly, opposite to the behavior found for an unmagnetized DEG, see. e.g., Refs. l3ll - [34ll . A decrease of 7£(i?) with 7 in 
the present case of a classical plasma is not attributable to the applied magnetic field because the field-free friction coefficient 
given by d40l shows a similar behavior (not shown in Fig. [3]). In a degenerate plasma an enhancement of the low-velocity SP 
with 7 is a quantum effect which is absent in our present study. For a DEG the domain of plasmon excitations is shifted towards 
smaller ion velocities dHHH]; this increases the SP in this velocity regime. But in the present case the domain of collective 
excitations is shifted towards higher velocities [47] and the friction coefficient decreases with 7. 

The resulting friction coefficient d47T > may be compared with Eq. (44) of Ref. Ill ill where the friction coefficient in the colli- 
sionless plasma contains an anomalous term ln(u e /u) vanishing at ■& — 0. The physical origin of such an anomalous friction 
coefficient may be traced to the spiral motion of the electrons along the magnetic field lines. These electrons naturally tend to 
couple strongly with long-wavelength fluctuations (i.e., small fcn) along the magnetic field. In addition, when such fluctuations 
are characterized by slow variation in time (i.e., small w = k • v), the contact time or the rate of energy exchange between the 
electrons and the fluctuations will be further enhanced. In a plasma, such low-frequency fluctuations are provided by a slow 
projectile ion. The above coupling can therefore be an efficient mechanism of energy exchange between the electrons and the 
projectile ion. At vanishing damping and in the limit of v — > 0, the frequency 10 = k • v —> tends to zero as well. The contact 
time ~ thus becomes infinite and the friction coefficient diverges. The collisions of the plasma particles play a stabilizing 
role since the fluctuations provided by the slow ion are damped. Thus at small velocities v — > the contact time is finite and is 
determined by ~ 7~ 1 . As a result Eq. ( |47| > does not contain a term like \n(v e /v) but behaves as 111(1/7) at vanishing damping. 



V. KINETIC VERSUS HYDRODYNAMIC APPROACH 



Using the theoretical results obtained in Sec. |IV] we present here some comparative analysis, looking for some contacts 
between our linear-response (kinetic) formulation and the previous hydrodynamic mode-coupling treatments based on the self- 
diffusion coefficients. As mentioned above such connection is established via Dufty-Berkovsky relation ( |49| >. The plasma is 
modeled as a collisional dielectric medium whose linear response function, within RTA, is given by Eqs. (Y2\ , <fl~8T> - (f20b with 
7 as a model damping parameter. In order to document the LIVSD physics highlighted by the relation ( |49| i, we first briefly 
pay attention to the unmagnetized B = limit. We consider it through the small e = l/(47rno e Ag) <C 1 plasma parameter 
approximation for the self-diffusion coefficient given by Sjogren et al. 114 811 . Employing Eqs. d35l l, (|4"TT i and d49b an inspection 
shows that at vanishing damping (7 = 0) the self-diffusion coefficient obtained from these formulas coincides with the result of 
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FIG. 4: (Color online) Proton LIVSD in a plasma with n e = 1.064 x 10 16 cm" 3 , T e = 1 eV (e = 0.02) in terms of Q e /w e . The lines with 
symbols in left and right panels represent parallel and transverse LIVSD, respectively. The solid lines were obtained from Eqs. d35t-J37b with 
# = (left) and ■& = tt/2 (right). 



Ref. |48] if the ionic charge number square Z 2 in Eq. ( f35T > is replaced by the quantity Z 2 — > P(Z), where 

, . / 1 \ 32Z 2 + 75V2Z + 50 

P{Z)=[Z+ — ) i-= . (50 

V V2J 104Z 2 + 111V2Z + 59 

For a proton projectile with Z = 1 this factor is P(l) = 1.003 and the agreement between both approaches is almost perfect. The 
factor d50l l which is nonlinear with respect to Z accounts for the nonlinear coupling between an incoming ion and the surrounding 
plasma B48I1 . However for highly charged ions with Z 3> 1 this factor increases linearly with Z, P(Z) = (4/13)Z, while more 
rigorous treatment shows that at strong ion-plasma coupling the energy loss of an ion scales with its charge approximately like 
Z r ' 0. 

Consider next the case of a magnetized and collisional plasma. For simplicity we consider electron-proton plasma and a proton 
as a projectile particle. Exploring first the moderately magnetized domain, fl e u e , one can explicit the field-free parallel and 
independent transverse diffusions 12411 . 

lf> = ^, Df = % (51) 

where v 2 — ksT e / m p , m p is the proton mass, and 7 C = w e e ln(l/ e) is the collision frequency in terms of the plasma parameter 
e, ri, — v p /n, e is the Larmor radius. Note that the collision frequency 7 C is related to the e-e collisional relaxation rate 7 ee as 
7 ee = ^/2/97T7 c (see Sec. [TlTb - The transverse diffusion coefficient given by Eq. (IBTt corresponds to a classical region, where 
Dj_ - B~ 2 , and is valid for 7c < il e < 0Aoj e Y(e) with Y(e) = [e 2 ^(l/e)]" 1 / 2 as explained in Ref. 0. 

With higher magnetic field values, f2 e /o/ e > 4Y(e), one reaches the transverse hydro-Bohm regime with D± ~ B^ 1 featuring 


D ± =D^ + J-e 2 [ln(l/e)] 3 / 2 . (52) 

To the intermediate plateau regime with D±_ ~ B° between transverse diffusion coefficients given by Eqs. ( Bil l and ( |52] >, 
corresponds the diffusion coefficient ll24ll valid at 0.4F(e) < Vl e /uj e < Y(e), 



D\ = D 



- n(o) 



i + 0.6c7 e a 



(53) 



with f = n e /to e . When electron diffusion is considered, the electron thermal velocity v e should be used in Eqs. d5Tl)-(|53"l) instead 
of v p . It is also important to stress that the quantitative predictions d5TT>-(|53T> of the mode coupling theory developed in Ref. ll24ll 
are strongly dependent on the values of the hydrodynamic cutoffs which, in contrast to the kinetic theory, are introduced linearly. 
A reliable estimate of the magnetic field and plasma parameter dependence of the cutoffs have been obtained, but their exact 
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values are not known. Consequently, the numerical coefficients in Eqs. (l5Tb-(T53T> are not precisely known. The exact coefficients 
can in principle be obtained by using kinetic theory. 

The friction coefficients S(v)/v (at v — > 0) calculated with the help of Eqs. (l49b and (BTl)-(l5"3l) are shown in Fig.[4]as the lines 
with symbols. In this figure the solid lines without symbols demonstrate the friction coefficients calculated from Eqs. d35T>-(f37T> 
with $ = (left panel) and $ = ir/2 (right panel) assuming, for consistancy, the same collision frequency 7 = 7 C as in Eqs. (T5TT >- 
653V As far as we can see, there are no fundamental contradictions between kinetic [Eqs. d35ll-(l3~7]>l and hydrodynamic [Eq. (|49| > 
and the first relation of Eq. (IBTl lI approaches for the parallel case, see Fig. |4] left panel. However, there are differences between 
the two approaches. It is well known that a kinetic equation approach contains more information about a physical system than 
a hydrodynamical approach. Indeed assuming a vanishing damping (7 = 0) and magnetic field (B = 0) from Eqs. d35l ) and 
( |4TT i at e <C 1 the ratio of the low-velocity SPs of both approaches is Sn ,kin/5'|| ,hyd — \/2. As discussed above the numerical 
coefficients in Eqs. (l5T1i-(l53l are not precisely known. Thus including the numerical factor y2 into denominator of the parallel 
diffusion coefficient in Eq. ( BTT l the agreement between both approaches becomes complete. Note that this is equivalent to the 
redefinition D\\ — »• Djj = Vp/j ee , where -f ee = y/2/9Trj c . 

An apparently large discrepancy is documented for the transverse situation (Fig. [4] right panel) where typically 
<S'_L,kin/<S'_L,hyd ~ [eln(l/e)] Q <C 1 with e <C 1 and a varies between 2 ^ a ^ 4 depending on the strength of B. The 
kinetic regime seems to be restricted to 0.08 = 7 c /w e < il e /uj e < 0AY(e) = 10.1 as explained in Ref. ll24ll . The discrepancy 
in the orthogonal case might be due to a different treatment of cutoffs in kinetic and hydrodynamic theories, i.e logarithmic vs 
linear. Actually, according to Ref. [24], the different hydro modes are normalized to distinct cutoffs. Upper hybrid ones are 
normalized to l/a e , a e being electron Larmor radius while low frequency modes are normalized to inverse mean free path \/t 
with I — v e /~f c . On the other hand in the extreme limit £l e S> co e < the only reasonable transverse cutoff should be l/a e which, 
for instance, in the kinetic treatment is included as ln(l/a e ). The basic physics involved in this orthogonal geometry pertains 
to kinetic theory when we rely on a collisional time while in the B — > oo limit leading to hydrodynamics, we incorporate the 
Larmor rotation of the charged particles, as well. 

VI. DISCUSSION AND OUTLOOK 

In this paper, we have presented a detailed investigation of the stopping power (SP) of low-velocity ions in a magnetized 
and collisional classical plasma. In the course of this study we have derived, among other things, some analytical results for the 
collision-inclusive linear response function for which the effect of collisions is taken into account in a number-conserving manner 
relaxation-time approximation (RTA) based on the Boltzmann-Poisson equations-the BGK approach. These analytical results at 
small projectile ion velocities go beyond those obtained in Refs. IU0l - [l5ll . After a quantitative introduction to the linear response 
function in Sec. [II] we have briefly studied the effect of magnetic field on relaxation rates due to collisions between different 
plasma species in Sec.[Ul]in which we have estimated contributions of close Coulomb collisions and long-range particle-wave 
interactions in the presence of a magnetic field on j a p. It has been shown that a magnetic field leads to an essential increase in 
collision rates. 

Theoretical calculations of SP based on the linear response theory within RTA are discussed in Sec. [IV] A number of limiting 
and asymptotic regimes of low-velocities and vanishing 7 have been studied. The theoretical expressions for SP derived in this 
section lead to a detailed presentation, in Sees. |IV| and[Vl of a collection of data through figures on SP of an ion. The results 
we have presented demonstrate that with regard to SP the difference between RTA and the usual linear response theory without 
damping is substantial. In particular, we have shown that the anomalous friction coefficient which behaves as S an (v)/v ~ 
ln(v e /v) at v <C v e obtained for a collisionless plasma ifTTll is now absent. Such a term arises due to an enhancement of energy 
exchange between plasma particles at v — > when the contact time ~ u^ 1 with u ~ k • v becomes very large. Thus collisions 
in a plasma play a stabilizing role for the energy exchange process and in the case of a collisional plasma the contact time is 
given by ~ 7 _1 . For low-velocity SP this yields S(v)/v ~ ln(l/7). Finally, in Sec.[IV]we have related the friction coefficient 
7£(#) at i? = to the parallel diffusion coefficient D\\ via Eq. ( |49l and at vanishing damping (7 — > 0), it has been shown that 
Du obtained from Eqs. (f35T> - (f37l> coincides with the result of Ref. 12511 . However, at finite 7 the equivalence of both approaches 
is violated due to damping. 

In Sec. [V] we have also compared the results obtained within our kinetic (Boltzmann-Poisson equations) approach with a sim- 
ple LIVSD expression ( l49| >, using transverse and parallel diffusion coefficients derived within a hydrodynamic mode-coupling 
theory 12411 . In the parallel case our results agree perfectly with this approach. There are discrepancies in the transverse situation 
and these have been discussed in Sec. [V] In the transverse case one may require an improved mode-coupling calculation. 

Going beyond the BGK approach which is based on the Boltzmann-Poisson equations we can envisage a number of avenues. 
These include (i) extending the number-conserving RTA to number-, momentum- and energy-conserving RTA. We note that this 
has been done previously in Ref. 14911 for a field-free case; (ii) studying SP in a degenerate electron gas (DEG) in the presence of 
a non-quantizing magnetic field, in a number-conserving RTA along the approach of Mermin [21] and Das B22I1 neither of whom 
considered a magnetic field; (iii) ion interaction with fluctuating plasma microfield and stochastic energy loss which may lead 
to a DB-like relation from first principles, and (iv) improved mode-coupling hydrodynamics which will result in an energy loss 
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calculation via DB relation and also exploring rigorous kinetic calculations of transport coefficients. Lastly, it may be mentioned 
that a quasiclassical model was studied by Das [50] which uses the Boltzmann-Poisson equations but also includes Landau 
quantization in a DEG. 
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Appendix A: Integral representation of the dielectric function 



In Sec. HI] we have derived the Bessel-function representation of the dielectric function. For some applications an integral 
representation of the dielectric function is desirable. For deriving the integral form of the dielectric function we rewrite the 
denominators of Eqs. Sl3[ and (fl4l using an integral 
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and summation formula J2^L-oo ^n( z ) eint — Jo(% z sin |) 12911 . Then the dielectric function may be alternatively represented 
in the form 
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where the argument of the Bessel functions is given by 
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In the above expressions by performing integration by parts, one finally obtains 
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In particular, for a Maxwellian isotropic distribution function (fTTl . Eqs. jA5i and (|A6t become 
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Here a e = v e /il e is the cyclotron radius of the electrons. The real and imaginary parts of the plasma dispersion function defined 
by Eq. ( fTSI ) can be found from Eq. jA7j and are given by 
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In this section we also derive an alternative integral representation for the collision frequency d27b in a magnetized plasma. 
Using an integral 
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and the summation formula (see, e.g., Ref. [29]) 
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from Eq. ( 1291 we obtain 



where 
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Let us recall that here a a — v a /tt a , v a , and fl a are the cyclotron radius, the thermal velocity and the cyclotron frequency of the 
plasma species a, respectively. 



Appendix B: Collision frequencies obtained from the Landau kinetic equation 



In this appendix we show that Eq. ( 1281 1 for the Coulomb logarithm is equivalent to the formula obtained by Silin in Ref. [46 
if one neglects the dynamical polarization of the plasma i.e. assuming e e i(k, u) = 1 in Eq. d28l . Note that the latter approach 
corresponds to the kinetic equation with the collisional integral in the form of Landau generalized for the case of a magnetized 
plasma. 

Using the integral representation ( IA14t of the function G a (k, u) we obtain 
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where U a (k, ut) is given by Eq. JA15t . Substituting this expression into Eq. (l28l for the Coulomb integral for the e-e collisions 
after k-integration we arrive at 
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where erf (x) is the error function. Similarly, for the e-i collisions we find 
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It is seen that Eq. ( IB2b with Eqs. ( IB3b - (IB5b and Eq. ( IB6b with dB7b coincide with the expressions derived in Ref. [46] on the 
basis of the Landau kinetic equation. 

In the case of vanishing or infinitely strong magnetic fields Eqs. (IB2b and iB6i become 



lnAo.ee = ln£ e , lnA , e j = lnf, 
lnAoo^e = iln^ e , In A^ = - ln£, 
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respectively. Thus, if the dynamical polarization effects are neglected, the Coulomb logarithm in the presence of a strong 
magnetic field is again given by | In Ao jeQ where, however, the quantity Ao ieQ is replaced by the usual one, i.e. r max /r m ; n . 
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